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Abstract. In the literature, the standard approach to finding bases of 
spaces of modular forms is via modular symbols and the homology of 
modular curves. By using the Eichler-Shimura isomorphism, a work by 
Wang shows how one can use a cohomological viewpoint to determine 
bases of spaces of cusp forms on To{N) of weight k > 2 and character \. 
It is interesting to look at the complexity of this alternative approach, 
and we do this for an explicit implementation of the algorithm suggested 
by Wang. 



When working in the field of modular forms it is often extremely useful 
to be able to work explicitly with spaces of modular forms, for instance in 
working with elliptic curves, testing conjectures, etc. 

To the author's knowledge, every available software package uses modular 
symbols to compute bases of spaces of modular forms, and in the literature 
this is the standard approach as well. A good reference for the computational 
aspects of this is Stein [7] . 

In this paper we look at another way of determining bases of spaces of cusp 
forms, by using a cohomological approach, based on the Eichler-Shimura iso- 
morphism, as done in Wang [8]. The theory behind this algorithm is worked 
through in [8], and what we do here is take an explicit implementation of 
this algorithm and analyze its complexity. 

Our main result is the following: 

Theorem 1. An upper bound on the theoretical complexity of determining 
a basis for Sk(To(N),x) via the cohomological approach described below is 



Finally we give two examples, where we work through the main steps 
of the algorithm. The first of these examples is the easier case of trivial 
character, while the second example with non-trivial character showcases 
some other aspects of the algorithm. 

1.1. Notation. Let A denote all 2 x 2 matrices with coefficients in Z, and 
let r = SL2(Z) be the matrices herein with determinant 1. We also define 
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0(N 3+e k 2+e (N + k A )) 



for e > 0. 
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We wish to determine a basis for the space of cusp forms of weight k > 2 
on Tq(N) with character x ( a Dirichlet character on (Z/iVZ)*). We denote 
by [i the index of Tq(N) in T. 

For J = e A , we define X (^o) = x(<0- 

As in [8], we let P^(Z/JVZ) denote 

{(c,d) G Z/iVZ x Z/AZ | gcd(c, d, N) = 1} 
modulo the relation 

(Ac, Ad) ~ xW(c, d), AG (Z/AZ)*, 

and we see that (Z/iVZ) is just P 1 (Z/AZ) if % is the trivial character. 
For 7 = ( a c b d ) G T we put 

7 = (c mod A, d mod A) G P^Z/iVZ), 

and it is easily checked that 707 = x(7o)7 for 70 G Tq(N) and 7 G T. 
We have an operation of A on 

' k-2 



3=0 



a 0i ■ ■ • j Ofc-2 G Z 



given by 



xV 5 = (ax + (ct + dyf' 1 ^, ( G A, 

and we get an action of Ao(A^) on M x (where the coefficient ring of M is 
extended with the values of x) by setting 5o- ?TJ = x(^o)(^o m )- 
We also look at the coinduced module W x of M x on T: 

W x = {w : T ->■ M x I 10(707) = 70-^(7) fo r 7o G r (A)}. 
We get an action of 5 G A on w G W x by setting 

r 5 ^M = / ' 7^A (A)r, 

l ' \5 Q .w{i), 7<5 = ^07', 5q G Ao(JV),7' G T, 

for 7 G r. 

For a matrix 5 G A, we denote by the submodule of W x invariant 
under the action of 5. 

We reserve /, S, Q and e for the following matrices: 

<=G !)•*=(; oV=(! "i i )"=(o 1 ? 

2. Algorithm and implementation 

2.1. Coset representatives. We start out by getting coset representatives 
for ro(iV) in T. Besides the A representatives ( ^ Tj 1 )>•••>( 1 n-1 ) > we a ^ so 
have (01) (if -/V is prime these are all of them) . 

As shown in Cremona [1], there is a bijection between coset representatives 
( a c b d ) of T (N) in T and elements (c, d) G P 1 (Z/iVZ). To get the remaining 
representatives we simply take the remaining elements (c,d) of P^Z/AZ) 
and lift these to a matrix ( ° b d ) in T via the Euclidean algorithm. 
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It is described in Section 2.2 of [1] how to efficiently determine P 1 (Z/A^Z). 
One starts out with the obvious elements (1, 0), (1, 0), . . . (1, N—l), and then 
look at elements (c, d), where c | N and d = 1, . . . , N — 1, and add it to the 
list if it is not equivalent to an element already on the list (one uses that 
two elements {c\,d\) and (02,^2) are equivalent if and only if c\d2 = C2d\ 
(mod N)). See also Section 8.7 of [7]. 

We will denote the coset representatives by 71, ... ,7^. 

2.2. Action of A on W x and relations matrix. From p. 103 of [8] we 

have the exact sequence (of complex vector spaces) 

— ► S k (T (N), x ) — H 1 (T,W X )+ — H\{T), W x ) + — 0, 
where T = ( q 1 ) • 

By p. 105 of [8] we have that H l (T, W x )+ S W x / (W£ +W$ +W£) , which 
is W x modulo the relations 

w + S.w = w + Q.w + Q 2 .w = w + e.w = 0. 

Therefore we need to be able determine a matrix representation of the 
action of a matrix in A on W x . 

An element of W x is determined by its values on the \x coset represen- 
tatives, and since M x is generated by the k — 1 homogenous monomials of 
degree k — 2, the space W x has fi{k — 1) generators. 

As a basis for W x we thus have the elements 

xiy k ~ 2 ~j , r = i, 




r j^i, 

with i = 1, . . . , /i and j = 0, . . . , k — 2. 

To determine the action of 5 € A on W x , we only need the action of 5.Wij 
on the coset representatives. If 7 r 5 ^ Ao(iV)r the action is 0, so we now 
assume that we can write ^y r 5 = <5o7 for some 5o G Ao(iV) and 7 € T. Since 
we have 7 = 70 7s for some 70 £ To(iV) and a coset representative 7 S , we 
replace 5q with <5q7o £ Aq(N), so that the action is given by 



5o.xiy h 2 J 'j s = i, 
0, s/i 



To get the action of <5 on t/;^ we need to run through the coset representa- 
tives j r , get the corresponding 5q € Aq(N) such that ^y r 5 = <5o7 s , and then 
compute the coefficients ao, . . . , a^-2 of the polynomial 

t=0 

for j = 0, . . . , /c — 2. These coefficients are then placed in the (s + j//)'th row 
and the (i + i//)'th columns (t = 0, . . . , fc — 2) of a /x(fe — 1) x fi(k — l)-matrix. 

This way we get matrix representations of the actions of I + S, I + Q + Q 2 
and I + e on W x , and + W% + W^) is then the nullspace of the 

resulting relations matrix (the relations matrix is the above three matrix 
representations stacked on top one another). 
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2.3. Dimension. To determine the dimension of the space of cusp forms, we 
need the dimension of if 1 ((T), W x ) + , that is the dimension of the subspace 
of -ff 1 ((T), W x ) invariant under the action of e. 

By Lemma 6 of [8] we have an isomorphism (of Q- vector spaces) 

H\{T),W X )= Q, 

s cusps 

and the action of e on ® s cusps Q is shown to be given by 

e -{ s ) = -x($o)W}, lj e = S oli, s = 7~ 1 oo, s' = l^oo, 
with i,j G {1, . . . 

Let us write e.{s} = c{s'}. Since e is an involution, we have £-{s'} = 
c {s}. If {s} = {V}, we have e.{s} = ±{s}, showing that {s} is in the cor- 
responding ±-space. If {s} ^ {s 1 }, we have {s} ±c{s'} in the corresponding 
±-space. 

Thus, if {s} / {■$'} we get elements of both ±-spaces, but when {s} = {s'} 
we get an element of only one of these spaces - the space corresponding to 
the sign of — x(^o)> i-e. we get at element of the +- space if and only if 
x{6o) = -1 when {s} = {s'}. 

We therefore need to determine when two cusps {7~ 1 oo} and { , y^ 1 oo} 
are equivalent, and this happens exactly when 7 m T u = 7o7n. for some 70 G 
Tq(N) and m£Z, This groups the coset representatives into classes (one 
for each cusp equivalence class). 

Thus e maps a cusp {s} to ±{s} (with s = ■j^ 1 oo) if and only if the coset 
representatives 7, and jj, satisfying jj£ = <5o7i, are in the same equivalence 
class, which is exactly when ^iT u = 707^ for some 70 G Tq(N) and mSZ. 

Since we only need one representative for each equivalence class, we choose 
for a given 7,, the unique representative 7^ satisfying jiT = 707,- for some 
70 G To(A r ), and write 7j£ = <$o7«) checking the sign of x(^o)- 

This way we find the part of the i-spaces coming from the case where 
e maps {s} to ±{s}. When this does not happen, we get elements of both 
±-spaces, and since the sum of the dimensions is we get the dimension 
of the +- space. 

The dimension of Sk(To(N), x) is the difference between the dimension of 
the nullspace of the relations matrix and the dimension just found. 

2.4. Hecke action and basis. Since Sk{To(N), x) is t ne kernel of the ho- 
momorphism 

H\T,W X ) + ^H\(T),W X ) + , 

we take elements in the kernel of this map and compute the corresponding 
(/-expansions until we have enough forms to generate Sk(To(N),x)- 

It is described in [8] how to choose elements in the kernel, and we will 
briefly recount this here. 

By Lemmas 3 and 4 in [8] we get that W x = M <g> (Z/iVZ) via the map 

w | ->$^7r 1 - TI, (7i)®7i- (*) 
i=i 
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Since H l ({T), W x ) = W x /(1 — T)W X , Lemma 6 of [8] gives an isomorphism 
(of Q-vector spaces) W x /(1 - T)W X = S cusps Q induced by the map M <g> 
pi(Z/iVZ)^0 scusps Q given by 

m® (c,d) h-> m(0, l){7 _1 .oo} = m(0,l){ - \ }, 7 = (c,d). 

Thus, the homomorphism H l {T,W x ) H X ({T),W X ) b ecomes a homo- 
morphism W X /(W X + W X ) — > S cus Q, and with the above isomorphisms 
this map is on W x = M ® P* (Z/iVZ) given by 

m ®( C) d)^m(0,l){-|}-m(l,0){§}. (f) 

For m <g> (c, d) to be in the kernel, we can use any (c, d) € P* (Z/iVZ) with 
m any (non-empty) linear combination of the monomials xy k ~ 3 , . . . , yx k ~ 3 if 
k > 2, and m = 1 if ft = 2 (in the weight-2 case we have to have x( c ) = x(d) 
as well). 

Since elements (c, d) € P^Z/iVZ) correspond bijectively to coset repre- 
sentatives (" ^) of ro(iV) in r, we can therefore represent elements of the 
kernel as m (8) 7r with r = 1, . . . , fx and m as above. 

We use the Heilbronn-Merel matrices (see Proposition 20 of Merel [5]) 



' c J.' 



ad — bc = n, a > b > 0, d > c > 



to determine the action of the Hecke operator T n on the elements of the 
kernel found above. 

The action of the Hecke operator T n on W X /(W^ + W x ) is given by 



T n x= ^ A 



x. 

AeH n 

Translating this through the isomorphisms above, the action on a kernel 
element m <8> 7r is 

T n (m <g>7 r )= x(M(^ m ) 8 Tr X i 

where we for each A £ -ff n write 7 r A = #o,A7r,4 with <5o,a £ Ao(iV) (terms 
where it is not possible to write ^ r A in this way are ignored). 

The coefficients of the polynomial x(<5o,A)(Am) are then saved in a vector 
t n , where the coefficient of x^y k ~ 2 ~^ is added to the r^ft — 1) — (ft — 2 — j)'th 
entry, as ^4 runs through i? n . 

If we wish to compute the basis up to exponent q M , we compute t n for 
n = 1, . . . , M, where we note that M has to be at least [y^rj j and we let t 
be the matrix whose n'th column is t n . 

We then multiply the nullspace matrix of the relations matrix found earlier 
with the matrix whose n'th column is t n , and denote the resulting matrix 
by B. If B has rank equal to the dimension of S/ t (To(N), x), we have found 
a basis (the leading rows of B). 

If B has rank less than the dimension, we choose another element m <8> 7 r 
in the kernel, compute the Hecke action ti,...,tjvf on this, multiply the 
nullspace matrix of the relations matrix with the resulting matrix, and get a 
matrix whose rows are concatenated to B, and we again compute the rank 
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of B. This procedure continues until we get as many linearly independent 
rows in B as the dimension of Sk(To(N),x)- 

Experimentation indicates that if we choose m = xy k ~ 3 + • • • + x k ~ 3 y (for 
k > 2) this procedure is likely to give a basis using just the first few 7 r 's. 

2.4.1. Determining the kernel. In the implementation described, we choose 
certain kernel elements and generate Fourier coefficients from these until we 
have enough forms to generate the space of cusp forms. 

However, there is no certainty that this will work, that is there is no 
guarantee that this approach will give enough linearly independent forms 
(even though the author has yet to see an example of this). 

Another approach (which is certain to work every time) is the following. 
From the nullspace of the relations matrix we have a basis for W X /(W^ + 

w x + W x)- B y usin S the isomorphism W x = M ® P* (Z/iVZ) given by (*), 
we can get the image of this basis on a quotient of MCEDP* (Z/iVZ), expressed 
in terms of the standard basis of M <8) P* (Z/iVZ). 

We use this basis to write up a matrix representation of the map (f ) on 
this quotient, and determine its nullspace (which has dimension equal to the 
dimension of Sk(To(N), x))- From the nullspace matrix we read off a basis, 
and we then write the kernel elements as linear combinations of the m®^ r . 
We now compute the Hecke action on these elements, which we know will 
generate enough forms to give a basis for Sk(To(N),x)- 

3. Complexity of implementation 

We always assume that the level is given via its prime factorization, i.e. 
no work is needed to find the divisors of N. 

We also use that a Dirichlet character on (Z/iVZ)* is defined via a lookup 
table, which takes O(N) to create, but we do not then need to worry about 
the cost of evaluating the character. 

3.1. Coset representatives. The number of coset representatives is 

/*=^II( 1 +i)- 

p\N 

Let n be the number of prime divisors of N, and let p±, . . . ,p n be the first 
n primes. By using Landau [3], p. 139, we find that 

n 

n( i +i)sn( i +^+-)=n( i -r I £ij>-ir I 

p\N p\N p\N i=l 

= 0(logp n ) = O(logn) = (log log TV), 

since N >2 n . We therefore have \i = 0(N log log N). 

We determine the coset representatives by looking at elements (c, d) , 
where c | N and d = 1, . . . , N — 1. By Theorem 315 of Hardy- Wright 
[2], the number of divisors of N is 0(N S ) for 5 > 0, and we therefore look 
at 0(N 1+S ) elements (c,d). 

Everytime we look at an element we check if it is equivalent to something 
already found, and this takes C(log 2 N) each time. Whenever we find a new 
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element we lift it to a matrix in T via the Euclidean algorithm, which also 
takes 0(log 2 N). This gives a total complexity of 0{N l+s log 2 N). 

3.2. Action of A on W x and relations matrix. All it takes to get -) r 5 on 
the form <5o7 s is to compute some greatest common divisors and run through 
the coset representatives to see which one works. All in all the complexity 
for this is 0(>log 2 N). 

The hardest part of computing the action on a matrix on a polynomial 
is the binomial coefficients that shows up when computing polynomial co- 
efficients. The cost of computing (™) is 0(m 2 log 2 n), and so a rough es- 
timate for computing the action of 5$ on the k — 1 monomials xiy k ~ 2 ~i 
is 0(k 4 log 2 k). This needs to be done for every j r , i.e. \i times, giving 
a total complexity of 0(fi 2 log 2 N + /ifc 4 log 2 k) for determining the matrix 
representation of 5 on W x . 

In can be noted that if we work over the finite field ¥ p instead of the 
integers, one can use a congruence first proved by Lucas [4] to obtain the 
complexity 0(p 2 log 2 p log A;) for the binomial coefficient computations. A 
more modern reference for the Lucas congruence is Stanley [6], p. 44. 

The relations matrix consists of the matrix representations of the actions 
of I + S, I + Q + Q 2 and I+s, and so is a matrix of size 3[i(k — 1) x (J,(k — 1). 
To compute the nullspace therefore takes 0(fi 3 k 3 ), and this is really the 
time-consuming function of this part. 

3.3. Dimension. We first build an array whose i'th entry is the index j 
of the coset representative satisfying jiT = 7o7j, as well as a similar array 
giving the index j of the coset representative satisfying 7,6 = <5o7j- Doing 
this takes C(/x 2 log 2 N). 

Next we determine to which cusp equivalence class each coset representa- 
tive belong, and computing along the way. The work needed is already 
done in the first array we created. 

We then choose a representative 7, of a cusp equivalence class and use the 
second array to find an equivalent representative jj satisfying jjE = <5o7i for 
some 5q G Ao(iV). We then compute x(#o) an d add 1 to the count of the 
corresponding ±- variable d±. 

The dimension of the +- space is then d + + (u^ — d + — d_)/2, and the 
dimension of Sk(Po(N),x) 1S the difference between the dimension of the 
nullspace of the relations matrix and the dimension of the +- space. 

The work in creating the arrays is by far the most work in this, so the 
complexity of this algorithm is 0(fi 2 log 2 N). 

3.4. Hecke action and basis. Even though we use the Heilbronn-Merel 
matrices in the implementation, we turn to [5] for the complexity analysis, 
since this paper gives another class of matrices which can be used instead 
of the Heilbronn-Merel matrices, and we have estimates on the size of these 
classes. 

In Section 3 of [5], Merel defines a set S n , where a matrix ("5) G A is in 
S n if it has determinant n and at least one of the following conditions are 
satisfied: 

• a > \b\, d > \c\, be > 0, 
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• b = 0, \c\ < d/2, 

• c = 0, \b\ < a/2. 

He also defines a set S' n , where \% b d ) £ A is in S' n if it has determinant n 
and one of the two following conditions are satisfied: 

• b = 0, \c\ = d/2, 

• c = 0, \b\ = a/2. 

It is easily seen that an upper bound for \S' n \ is 2a(n), where c(n) is the 
sum of the positive divisors of n, and from p. 85 of [5] we have, as n — > oo, 

K I 121 °g 2 ( \ 1 

o n ~ 5 — cr(n) log n. 

By Theorem 322 of [2] we have cr(n) = 0(n 1+s ) for S > 0, so that we have 

\S n US' n \ = 0(n 1+s logn) = 0(n 1+e ), 

for e > 0, since logn = 0(n s ) for 5' > 0. 

We need to compute S n and S' n for n = 1, . . . , |_ J ( or rnore n if we want 
higher precision), so 0(fi 2Jre k 2+e ) matrices are needed to compute the Hecke 
action. 

We want to compute the action of T n on m®^ r , for m a linear combination 
of the monomials xy k ~ 3 , . . . , x k ~ 3 y if k > 2 and m = 1 if k = 2. To do this 
we write, for each i £ 5 n U5 n , 7 r A = (Jo^Tr-^, with 5o t A £ Ao(iV), and this 
can be done in C(//log 2 iV) for each A. We then have 

T n {m®~i r )= ^2 x{So,A)( A - m ) ® 7r A + ^ X] x(^o,a)(-4-"i) ® 7r A5 

and as mentioned earlier, the complexity of determining the action of A on 
a linear combination of all possible monomials is 0(/c 4 log 2 k). 

The coefficients of x{3o,A)(A.m) are added for each A and saved in a 
vector t n , with indices depending on the monomial xiy k ~ 2 ~i and index r^. 
Determining all necessary t n are done in 0(fi 2+e k 2+e (fi log 2 iV + A; 4 log 2 k)). 

Multiplying the nullspace matrix of the relations matrix with the matrix 
whose n'th column is t n , takes C(/x 3 /c 3 ), which is less than the complexity 
of determining the Hecke action. 

The resulting matrix is the basis matrix if it has rank equal to the dimen- 
sion of Sk(TQ(N),x)- We therefore do Gaussian elimination and compute 
the rank to see if we are done. If not, we choose another coset representative 
7 r , get the resulting t n 's of the Hecke action on m<g>7r, multiply the nullspace 
matrix of the relations matrix with the resulting matrix, and get a matrix 
whose rows are concatenated to B. We again do Gaussian elimination and 
compute the rank, and this is repeated until we get the right rank. If we 
run through all the coset representatives without getting the right rank, we 
can try with another m. 

Gaussian elimination is done in 0(fi 3 k 3 ), and is therefore insignificant 
compared to the computation of the Hecke action. 

For a given m, this procedure is repeated at most ji times, but exper- 
imentation indicates that it is likely to finish much sooner if m is chosen 
to be the sum of all possible monomials (in the case of k > 2). We see 
that the computation of the Hecke action is by far the hardest part of this 
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basis determination, and we therefore get a total theoretical complexity of 
0{^ +e k 2+e {^i\og 2 N + k 4 log 2 k)). 

Since \i = 0(N log log N), and hence \f = C(iV e ), we find that 

0(fi 3+e k 2+€ (filog 2 N + k 4 log 2 k)) = 0(n A+e k 2+e log 2 N + fi 3+e k 6+e log 2 k) 

= 0( f i 4+ ^'k 2+t + ^k 6+ ^') 

= 0(N 3+e+e 'k 2+e+e '(N + k A )), 

using that log 2 a = 0(a e ') for e' > (and a either N or k), and this is 
Theorem 1 after replacing e + e' with e. 

4. Examples 

We use this implementation to determine bases for two spaces of cusp 
forms, one with trivial character and one with non-trivial character. 

The first example gives more detail, while the second highlights an aspect 
which only happens in the case of non-trivial character. 

4.1. S4(To(25)). We start out by getting the coset representatives, and we 
do this by determining P 1 (Z/25Z) in the way we described earlier. Thus we 
get the 30 elements 

(1,0),..., (1,24), (0,1), (5,1), (5, 2), (5, 3), (5, 4), 

and these are lifted to matrices in T via the Euclidean algorithm: 

-l\ fo -l\ fl 

1 J ' " ' ' \1 24 J ' ^0 1 

1 °) (~ 2 ( 2 (~ l 

5 y'V 5 2 /v 5 3 /v 5 4 

We denote the representatives by 71, . . . ,730, e.g. 726 = I- 

Next we determine the nullspace of the relations matrix. Since the matrix 
representations of I + S, I + Q + Q 2 and I + e all are 90 x 90-matrices, we 
do not write these up here. After bringing the matrix on echelon form and 
deleting zero rows, the nullspace matrix of the relations matrix is a 7 x 90- 
matrix. 

Next we compute the dimension of the +-space. We therefore build an 
array whose i'th entry is the index j of the coset representative jj satisfying 
7iT = 7o7j for some 70 € IV This becomes 

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 



J 



6 7 8 9 10 11 12 13 14 15 16 



i 


16 


17 


18 


19 


20 


21 


22 


23 


24 


25 


26 


27 


28 


29 


30 


j 


17 


18 


19 


20 


21 


22 


23 


24 


25 


1 


26 


27 


28 


29 


30 



This array shows that 71, ... , 725 are in the same cusp equivalence class, 
and the rest are in their own class. All in all, we see that To(25) has 6 cusps, 
represented by 71, 7 26 , 727, 728 , 729 and 730 . 

We also build an array whose i'th entry is the index j of the coset repre- 
sentative 7j satisfying 7^6 = <5o7i for some 5q G Aq(N). This becomes: 
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i 


1 


2 


3 


4 


5 


6 


7 


8 9 


10 


11 


12 


13 


14 


15 


3 


1 


25 


24 


23 


22 


21 


20 


19 18 


17 


16 


15 


14 


13 


12 


i 


16 


17 


18 


19 


20 


21 


22 


23 24 


25 


26 


27 


28 


29 


30 


3 


11 


10 


9 


8 


7 


6 


5 


4 3 


2 


26 


30 


29 


28 


27 



We now take each of the six representatives 7$, match them with the 
corresponding jj from this table (checking that 7j is in the same equivalence 
class). We do not need to compute the <5o's since we have trivial character, 
and therefore always will get elements of the — space. 

From the table we see that only 71 and 726 give rise to jj's in the same 
cusp equivalence class, and we therefore get 2 elements of the — space. 

Since there are 6 cusps, this means that the remaining 4 are split evenly 
between the +- and — spaces, giving that the dimension of the +-space is 
2. 

From this we get that the dimension of S^(T q(25)) is 7 — 2 = 5 (nullspace 
dimension minus +-space dimension). 

Finally we compute the Hecke action t±, . . . , tM for M = L^rJ = 10. 

Each coset representative 7 r gives rise to a kernel element xy ® 7 r , and 
to compute t n for xy ® 7 r we compute A.xy for all A £ H n , keeping track 
of the index r& of ^ r A = <5o,A7A r - 

In the case of n = 3 we have 

Hz = { (0 3) ' (l 3) ' (2 3) ' (l 2) ' (0 l) ' (0 l) ' (0 l) } ' 

and we denote these A\, . . . , Ay. 

We need to compute up to r = 4 since the forms generated by using r < 3 
only spans a 4-dimensional space. As a matter of fact, the forms generated 
by using just xy <g> 74 span the whole space, and we write a table of indices 
Ti with respect to 74, as well as the action of Ai on xy: 



i 


1 


2 


3 


4 


n 


10 


22 


13 


28 


Ai.xy 


3xy 


?>xy + x 2 


3xy + 2x 2 


2y 2 + 5xy + 2x 2 


i 


5 


6 


7 




n 


2 


19 


11 




Ai.xy 


3xy 


y z + 3xy 


2y 2 + 3xy 





Thus we get that 
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T 3 (xy (8)74) = y~] Aj.xy ® 7 ri , 
i=l 

and we put in the coefficients in a vector £3 (where (t3) s r r ._u + j is the coef- 
ficient of x- ? y 2 ~ J in Ai.xy). 

We do this for all n = 1, . . . , 10, and build a matrix with the t n as columns. 
We then multiply the nullspace matrix with this matrix, and get (after 
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removing zero rows and putting the matrix on echelon form) 



/I 























1 


o\ 





1 














-1 


-1 














1 











1 


-2 

















1 





-1 








-3 





v° 











1 














-v 



which has rank 5, and the rows therefore form a basis for S4(Tq(25)). Thus 
the standard basis (up to g 10 ) of this space is 

q + q 9 , q 2 - q 7 - q 8 , q 3 + q 7 - 2q S , <? 4 - q 6 - 3q 9 , q 5 - 4q 10 . 

4.2. 5s(ro(12), (jtj)). We start out by getting the coset representatives, and 
we again do this by determining P 1 (Z/12Z). We get the 24 elements 

(1,0),..., (1,11), (0,1), (2,1), (2, 3), (2, 5), (3,1), 
(3, 2), (3, 4), (3, 7), (4,1), (4, 3), (4, 5), (6,1), 

and these are lifted to matrices in T via the Euclidean algorithm: 

(0 -l\ (0 -1\ (1 0\ (1 0\ (1 l\ (I 2\ (1 0\ 

\i o y ' • • • ' Vi ny'v V'V 2 v V 2 3 /'V 2 5 /'V 3 V' 
/-l -l\ A A /l 2\ A 0\ /-I -A A A A 0\ 
^3 2j'^3 V'V 3 V V 4 V V 4 3 /V 4 5 /V 6 V ' 

We denote the representatives by 71, . . . ,724- 

We again determine the relations matrix (of size 184 x 96) and its null- 
space, which in this case has dimension 8. 

Just as in the last example we choose a coset representative from each 
cusp equivalence class, and we get 6 cusps represented by 71, 713, 714, 717, 
721 and 724- 

We now take each of these 7$, match them with the corresponding 7^ 
satisfying ■jjE = <5o7i (checking that jj is in the same equivalence class as ji, 
which they all are in this case), and compute x(<5o) to see to which of the 
±-spaces they correspond. We summarize the results in the follwing table: 



i 


1 


13 


14 


17 


21 


24 


3 


1 


13 


16 


20 


23 


24 


So 






( -5 2\ 
1-12 5 ) 


( -7 2\ 
\ -24 7 J 


( -24 5 ) 


U?) 


x(6o) 


-1 


1 


-1 


1 


-1 


1 


Space 


+ 




+ 




+ 





From this we get that the dimension of 55(^(12), (yg)) is 8 — 3 = 5 
(nullspace dimension minus +- space dimension). 

As before we compute the Hecke action ti, . . . ,ijvr f° r M = L^LT-I = ^> 
and in this case it is enough to use the kernel element (xy 2 + xy 2 ) (g) 71 to 
generate a basis. 
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The matrix we get after multiplying the nullspace matrix with the matrix 
having the t n as columns (and removing zero rows and putting it on echelon 
form) is 



fl 

















-4 





-27 


o\ 





1 











-3 





-8 














1 











-10 





12 














1 





-3 











6 













1 





-5 





9 


V 



We thus get that the standard basis (up to q 10 ) of S5(Tq(12), (j^)) is 
q-4q 7 -27q 9 , q 2 -3q 6 -8q 8 , q 3 - I0q 7 + 12q 9 , q A -3q 6 + 6q 10 , q 5 -5q 7 + 9q 9 . 
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